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We report on specific heat measurements on clean overdoped BaFe 2 (Asi-xPx )2 single crystals performed 
with a high resolution membrane-based nanocalorimeter. A nonzero residual electronic specific heat coefficient 
at zero temperature 7 r = C /T\t^o is seen for all doping compositions, indicating a considerable fraction of the 
Fermi surface ungapped or having very deep minima. The remaining superconducting electronic specific heat is 
analyzed through a two-band s-wave a model in order to investigate the gap structure. Close to optimal doping 
we detect a single zero-temperature gap of Ao ~ 5.3 meV, corresponding to Ao/ksTc ~ 2.2. Increasing 
the phosphorus concentration x, the main gap reduces till a value of Aq ~ 1.9 meV for x = 0.55 and a 
second weaker gap becomes evident. From the magnetic field effect on 7 r, all samples however show similar 
behavior [7r(F7) — 7r(7T = 0) oc Ff", with n between 0.6 and 0.7]. This indicates that, despite a considerable 
redistribution of the gap weights, the total degree of gap anisotropy does not change drastically with doping. 

PACS numbers: 74.25.Bt, 74.70.Xa 


I. INTRODUCTION 

The symmetry and stiTicture of the superconducting gM 
in iron-based superconductors is a highly debated topioM 
In contrast with cuprates, where a tf-wave symmetry is 
predominanPS, iron-based superconductors present a rela¬ 
tively wide range of possible sc enari os. The leading hypothe¬ 
sis is a multi-band s-t symmetr^DHS where the order parame¬ 
ter assumes opposite sign on different sheets of the Fermi sur¬ 
face, but remains relatively constant in amplitude along them 
(no nodes). Despite the fact that many iron-based supercon¬ 
ductors present a gap structure compatible with a nodeless 
nodal behavior has been observed in several compounds^G^l 
Accidental nodes, not due to the gap symmetry but to strong 
variations of the gap amplitude along a Fermi surface sheet, 
have then been taken into account in theoretical models to 
reconcile the app arent symmetry with zeroes in the su¬ 
perconducting gapl®3. A clear picture has however not been 
achieved yet and more experimental and theoretical efforts are 
required. The isovalently doped system BaFe2(Asi_a;P3;)2 of 
the 122 family is particularly interesting as signs of nodal be¬ 
havior have been detected even at optimal dopin^MUl This 
is in stark contrast with its hole- and electron-doped coun¬ 
terparts which, despite showing similar phase diagrams and 
critical temperatures, are believed to be nodeless at least at 
optimal dopings®. Specific heat measurem ents hav e been 
performed on the hole-doped Bai_a;Ka;Fe2AsJ^ ^I^^I^^ I and on 
the electron-doped Ba(Fei_a;Co3^)2As5p3^. While some fea¬ 
tures are common to both compounds, e.g., multi-gap behav¬ 
ior and strong-coupling values of the main gap amplitude, oth¬ 
ers are substantially different, such as the values of the resid¬ 
ual specific heat coefficient 7^ and the magnetic field- and 
doping dependence of the specific heat. 

In this work, we study the specific heat of 
BaFe2(Asi_a;Pa;)2 single crystals in the overdoped regime. 
The low temperature electronic specific heat provides infor¬ 
mation about the gap amplitudes, while its magnetic field 


dependence at low temperatures reveals the gap anisotropy. 
A sizable value of the residual electronic specific heat 
coefficient 7r, amounting to 17% near optimum doping, 
growing to ~ 30% with doping, shows that a considerable 
part of the Fermi surface is ungapped or presents broad deep 
minima. The remaining superconducting specific heat is 
fitted to a two-band a-model, which is found to represent the 
experimental data well. Close to optimal doping, the best fit 
is obtained with a single gap function with a zero-temperature 
energy gap Ag ~ 5.3 meV. However, a second gap becomes 
evident as x increases. The magnetic field dependence of the 
zero temperature specific heat coefficient reveals a sublinear 
behavior [A7 = 7r(iT) — 7r(7F = 0) cx iT(0'6-0-7)j^ largely 
independent of doping. 


II. EXPERIMENTAL DETAILS 

The high-purity BaFe2(Asi_xPx)2 single crystals inves¬ 
tigated in this work were grown with a self-flux method^. 
Three crystals with composition x = 0.32, 0.50 and 0.55 were 
selected and cleaved in order to obtain plate-like samples of 
side ^ 100 — 200 pm. Optical microscopy observation corn- 
firmed that all samples had shiny surfaces with no identifi¬ 
able secondary phase inclusions, x = 0.32 cotTesponds to 
near optimal doping with = 28.4 K, while x = 0.50 and 
x = 0.55 are in the overdoped regime with = 18.2 K and 
12.5 K, respectively. Specific heat was measured with a dif¬ 
ferential membrane-based nanocalorimeter applyi ng an AC- 
method with phase stabilized frequency feedbaclPSESl The 
sample side cell of a typical calorimeter device is illustrated 
in Fig. [TJa). The active area is a stack of thin Aims in the 
center of a Si3N4 membrane, whose key elements are a GeAu 
thermometer and Ti (AC and DC) heaters. A minute amount 
of Apiezon-N grease was used to attach the crystals to the 
calorimeter, as shown in Fig. [^b). The grease specific heat 
was measured separately in order to subsequently extract the 
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FIG. 1. (a) Optical microscope image of the sample side of a typical 
nanocalorimeter, (b) SEM image of the sample cell, with the most 
overdoped sample (Tc = 12.5 K) attached to it. (c) Specific heat 
plotted as C/T of the Tc = 28.4 K sample as a function of temper¬ 
ature T. Red circles represent the experimental data, while the black 
curve is an interpolation of the normal state signal with a Debye- 
Sommerfeld function, (d) Temperature dependence of the electronic 
specific heat Ce at low temperatures shown as Ce/T for the same 
sample as in (c). 


intrinsic sample response. 


III. RESULTS AND DISCUSSION 


The temperature dependence of the specific heat plotted as 
C/T of the Tc = 28.4 K crystal is reported in Fig.[^c). The 
transition to the superconducting state is apparent as a small 
peak at Tc. The specific heat above the transition is fitted with 
a Debye-Sommerfeld function and extended below the transi¬ 
tion with the requirement of entropy conservation at Tc. The 
electronic specific heat is obtained by subtracting the phonon 
contribution from the total specific heat and is shown as Ce/T 
in Fig.[TJd) for the same sample. A similar procedure was ap¬ 
plied for the Tc = 18.2 K and 12.5 K samples. In the latter 
case, with a relatively low upper critical field 11 ^ 2 , it is possi¬ 
ble to verify that the normal state when applying a magnetic 
field of 5 T actually corresponds to the one calculated with the 
Debye-Sommerfeld function. 

In all studied samples, a sizable residual specific heat 71 
is found, ranging from 17% of the normal state coefficient 
7 n for Tc = 28.4 K to 31 and 28% for the Tc = 18.2 K 
and Tc = 12.5 K samples, respectively. This term is due 
to the presence of non-superconducting quasiparticles. Since 
the crystals are nicely shaped with freshly cleaved surfaces, 
macroscopic secondary phases are unlikely to be the cause of 
this high 7 r. Therefore, we believe that the residual term is 


due to a part of the Fermi surface being ungapped or present¬ 
ing deep broad minima in the order parameter suppressed by 
weak disorder. Values of 7 r/ 7 n in the range of 10 — 20% 
have been reported for Co-doped samples close to optimal 
dopin^21l25| ( 7 ^/ 7 ^ ^ 5% at optimal doping), much higher 
than for K- dope d samples, where 7 r/ 7 n accounts to only a 
few percenPmU. This is in agreement with a fully gapped 
state detected in K-doped samples over a wide range of the 
phase diagratrl^, while nodal behavior is detected from ther¬ 
mal conductivity measurements in Co-doped samples as soon 
as X moves away from optimal doping. The rather high value 
of 7 r/ 7 n in our samples is in agreement with nodes/minima 
measure d fro m angle-resolved ^otoemission spectros copy 
(ARPESj®?, penetration deptlP^, thermal conductivitjWll^ 
and nuclear magnetic resonanc^i^. The absolute value of 7 ^ 
in our samples does not show an increase with doping as re¬ 
ported for Co-doped Ba-12i231. Since quantum oscillatio ns are 
observed only in overdoped BaFe 2 (Asi_xPx )2 sample ^221311 
it is likely that phosphorus doping decreases the disorder in 
the crystal, avoiding a strong increase in 7 r. 

In order to obtain the superconducting contribution Ces/T 
the residual term 71 . is subtracted from the total C/T. The re¬ 
sult is shown in Fig. for the three doping levels analyzed. 
The data are presented as a function of T/Tc and normal¬ 
ized by ( 7 n — 71 ) in order to allow a direct comparison be- 



FIG. 2. Temperature dependence of the reduced electronic specific 
heat Ce/T, subtracted by 7^ and normalized by (7n — 7 r). Data 
points are represented by colored dots, the fits according to a two- 
gap a-model by black lines. The fitting parameters are reported next 
to each curve, (a) T = 28.4 K. (b) T = 18.2 K. (c) T = 12.5 K. 
(d) Resume of the 3 data sets. 
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TABLE I. Superconducting parameters for BaFe2(Asi_xPx)2 ob¬ 
tained from electronic specific heat. 


Property 

Unit 


TcfK) 


28.4 

18.2 

12.5 

7 

n 

mj/molK^ 

28.0 

20.2 

19.2 

7r/7n 

% 

17 

31 

28 

ai 


2.15 

2.22 

1.8 

012 


- 

0.85 

0.52 

Ai 

meV 

5.3 

3.5 

1.9 

A2 

meV 

- 

1.3 

0.6 


n 

% 

100 

93 

71 

n 

^avg 

meV 

5.3 

3.4 

1.7 

r 

^cond 

meV 

4.6 

2.9 

1.7 

a 

% 


0.64(6) 

0.66(1) 

0.68(2) 


“ Weight associated with the first gap 
** Average gap obtained according to Eq. ( ^ 

Average gap obtained according to Eq. ( 

Erom fittings of A'y{H) to the function*A 7 = A ■ (fioH)" 


tween them. Cgs/T is then fitted to a two-b and implemen¬ 
tation of the phenomenological a- model 22 l 221 

which assumes 

full s-wave gaps both closing at Tc and a BCS temperature 
dependence. ARPES measureme nts in dicate that this is in¬ 
deed a reasonably good assumptiorP^Ii^ and as a consequence 
the a-model is able to give a good description of the specific 
heat data. An alternative approach would be the 7 -model of 
Kogan, Martin and Prozorov23 which can be used for a gen¬ 
eral temperature dependence of the gap, but which requires 
knowledge about the Fermi velocity and density of states on 
the different bands. The electronic specific heat according to 
the ct-model is calculated a^^ 


Cesjt) 
(7n -lr)Tc 


6 a^ 

TT^f 


nOO 

/ /(I-/) 

Jo 



ldA2\ 

2~dr j 


de, 

( 1 ) 


where a = Ao/fcBTc is an adjustable parameter, t = T/T^ 
e = e/Ao is the normalized single-particle energy, E = 
\/? + A^ is the normalized energy of elementary quasipar¬ 


ticle excitations, f{E,t,a) = ex-p{aE/t) + l 


1 -1 


is the 


Fermi-Dirac distribution function, and A(f) = A(f)/Ao is 
the order parameter normalized by itz zero temperature value. 
The total electronic specific heat Ces is then considered as the 
sum of two independent contributions Ci and C 2 , given by 
two different a values, ai and a 2 - A weight r/i is associated 
with Cl and r /2 = 1 — tyi with C 2 . A routine was imple¬ 
mented in order to vary the three free parameters ai, 02 and 
r]i to minimize the root-mean-square deviation between the 
data and the model function. 

The best fitting parameters obtained are presented in 
Fig.|^a)-(c) for each curve and summarized in Table Their 
uncertainties are estimated to ^ 5%. All fitting curves allow 
a fairly good representation of the data. For — 28.4 K 
the experimental curve clearly saturates at low temperatures 
{T/Tc < 0.25), leading to the conclusion that no minor gap 
should contribute, at least within our measurement resolution. 
The best fit is in fact given by a single gap function, with a = 


2.15 (A = 5.3 meV). For Tc = 18.2 K, the curve no longer 
saturates at low temperatures and a contribution from a second 
smaller gap has to be considered. The fitting routine gives a 
main gap ai = 2.22 (Ai = 3.5 meV), with a considerable 
weight rji = 93% and a minor a 2 = 0.85 (A 2 = 1.3 meV) 
for the remaining 772 = 7%. For Tc = 12.5 K, the main peak 
reduces considerably in amplitude (ai = 1.8, with an asso¬ 
ciated Ai = 1.9 meV). Moreover, the second gap contribu¬ 
tion is much more pronounced (772 = 29%). Its amplitude is 
A 2 = 0.6meV (a 2 = 0.52). The ais associated with the 
main gap are all higher than the BCS value aecs = 1-764, 
indicating strong-coupling. However, the Tc — 12.5 K sam¬ 
ple presents an ai much closer to aecs than the other two 
samples. A weighted average of the gap amplitudes Aavg is 
calculated from the a-model results for the three doping com¬ 
positions: 

Aavg = ^Ai^ • 771 + A2^ • 772 ( 2 ) 

The results are listed in Table The average value is dom¬ 
inated by the main gap in all samples. The obtained AavgS 
are compared to the gap amplitudes calculated from the BCS- 
style relation for the condensation energy 

AF = (3) 

7n 

where N{Ef) is the density of states at the Fermi energy, 
N{Ef) = 37n/7r^fc0. AF is obtained by integration of the 
AC (T) curve between 0 and tJ^, where AC is the difference 
in specific heat between the superconducting and the normal 
state. The gap values obtained from AF are reported in Ta¬ 
ble [I] as well. The most overdoped sample, which has an en¬ 
ergy gap ratio very close to the BCS value, shows very good 
agreement between Aavg and Acond- For the other two sam¬ 
ples, the average gaps are around 15% higher than Acond- 
This behavior is likely due to strong-coupling at these dop¬ 
ings, for which Eq. 0, valid under the assumption of weak- 
coupling, underestimates A. It is in fact theoretically ex¬ 
pected that strong-coupling superconductors show a conden¬ 
sation energy lower than that expected in the BCS cas d^^F? ! 



FIG. 3. Field dependence of A7 = 7r(77) — 7r(77 = 0) for Tc = 
28.4 K (red squares), E = 18.2 K (blue triangles) and E = 12.5 K 
(green circles). The curves are fits of the type A7 = A - (po77)"> 
with n given in Tablej^for the different samples. 
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A lower AF then corresponds to a lower Acond- The a- 
model instead phenomenologically takes into account strong¬ 
coupling effects and gives a better representation of the exper¬ 
imental data. 

The gap amplitude close to optimal doping, Aavg = 
5.3 meV, can be compared with the ARPES values of Refs. 
iMl and ca, which report ^ 5 — 8 meV for the gap on 
the hole pockets and ~ 7 — 8 meV on the electron pockets. 
The values seem in good agreement considering that we are 
slightly out of optimally doped conditions and that the spe¬ 
cific heat signal is more sensitive to the hole pockets (see OSl 
and references therein). P overdoping has the effect of reduc¬ 
ing the main superconducting gap and increasing the weight 
of the second gap. In order to better visualize the evolution of 
the gap with doping, experimental data for the three compo¬ 
sitions studied are plotted together in Fig. |^d). The behavior 
is different in comparison with Co-doped sampleJ^Hll where 
the two gaps have fairly constant weights throughout the en¬ 
tire doping range. It is instead in agreement with K-doped 
sampIegSS where the weight of the second gap is increasing 
with doping, as in the present case. 

The field dependence of the zero temperature specific heat 
coefficient 7 r is analyzed in order to extract information on the 
gap anisotropy. Specific heat was measured in magnetic fields 
up to 5 T for all samples. 5 T corresponds to 11%, 39%, and 
93% of 7 Jc 2 for the samples with Tc = 28.4 K, 18.2 K, and 
12.5 K, respectively. The zero temperature 7 ^ values are ob¬ 
tained by linear extrapolations of the low-temperature Ce/T 
curves. The resulting 7 r(Tf) values are shown in Fig. [^as 
a function of applied field figH for the three doping com¬ 
positions. The curves are vertically shifted down to zero by 
subtracting the zero field 7 r(i 7 = 0 ) from all 7 r(Tf) values. 
A nodeless s-wave order parameter is expected to give a lin¬ 
ear dependence of 7 ^ as a function of fiel(fl3, while a d-wave 
order parameter a square root dependence with fielcPI. An 
intermediate behavior is generally interpreted in iron-based 
superconducto rs as due to the presence of gaps with differ¬ 
ent amplitudeJ^^Illl. The gap imbalance becomes more pro¬ 
nounced as 7r(i7) deviates further from a linear relation. 
All curves were htted with a function of the type A 7 = 
A ■ (/xoiT)", with A and n being the htting parameters. The 
obtained exponents n are reported in Table |I] The exponent 
is in the range of 0.6 to 0.7 for all three samples. Previous 
measurements at optimal dopin^Sl show as well a very sim¬ 
ilar field dependence for fields up to 4T, while a crossover 
to a linear behavior is found at higher fields. This behavior 


is interpreted in terms of a double gap system, in which the 
Volovik-like trend at low fields is due to a strongly anisotropic 
gap, while the linear component at high fields is due to a sec¬ 
ond isotropic gap. Even if we see no sign of a second gap 
from the temperature dependence of the specific heat of the 
Tc = 28.4 K sample as shown in Fig.|^a), such a gap should 
be expected in order to explain the held dependence of A 7 . 
Its weight is however at most about 5%, given by the uncer¬ 
tainty of the htting parameters. The exponent n displays no 
clear trend with doping. This shows that, despite different 
weights of the gaps along the doping range, the total degree of 
anisotropy, as measured by the held dependence of A 7 , stays 
constant. This is in qualitative agreement with the theory of 
BanJ^, which shows that for a two-band s± state in the pres¬ 
ence of impurity scattering the held dependence of A 7 mainly 
depends on the ratio between the two gap amplitudes A 2 /A 1 
and not on their weights. 

IV. CONCLUSIONS 

In conclusion, from high-resolution specihc heat measure¬ 
ments on overdoped BaFe 2 (Asi_xPx )2 single-crystals we 
observe: (i) A sizable residual term j,- at all dopings, sign of a 
considerable part of the Fermi surface presenting no supercon¬ 
ducting gap or very deep minima, (ii) A main gap amplitude 
in the strong-coupling limit (a = Aq/ZcbTc between 1.8 and 
2.22), in agreement with ARPES measurements near optimal 
doping, (iii) A reduction of the main gap weight with doping, 
which passes from dominating the specihc heat signal near 
optimal doping to a 71% of the total weight for Tc = 12.5 K. 
(iv) A high gap anisotropy, constant at all dopings. 

Note: During preparation of this paper, we became aware of 
a recent specihc heat studjHS which shows that a single-band 
anisotropic gap ht gives a good representation of the experi¬ 
mental data as well. 
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